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The optical spectra in the family of multiferroic manganites _RMn03 is a great puzzle. Current 
models can not explain the fact that two strong electromagnons are present in the non-collinear 
spin cycloidal phase, with only one electromagnon surviving the transition into the collinear spin 
sinusoidal phase. We show that this is a signature of the presence of anomalous magnetoelectric 
coupling that breaks rotational invariance in spin space and generates oscillatory polarization in the 
ground state. 
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I. INTRODUCTION 

In multiferroic materials magnetic and electric orders 
coexist simultaneously and the coupling between spin 
and charge degrees of freedom gives rise to a wide range 
of magnetoelectric phenomen£p"2l Recent research has 
centered on the origin and symmetry of magnetoelectric 
coupling. The crucial question is how the coupling be- 
tween two spins depends on electric field: 

Hxae = ^ ] j^7im(E)S n • S m + D nm (E) • S rl X S m 

nm 

(E)-S m . (1) 

Here S„ and S m are spins at lattice sites R ra and R m , 
and the electric field E can be either internal, i.e., from 
the electric polarization in the material, or external, as 
is the case of incident light. The first two interactions 
in the right hand side of Eq. are well understood. 
The first one, exchange interaction J, is electric-field de- 
pendent because atomic positions are modulated by E 
(the phenomena of magnetostriction). The second one, 
the Dzyaloshinskii-Moriya (DM) interaction, is first or- 
der in spin-orbit coupling and is antisymmetric under 
spin interchange. The third and final interaction, the 
anomalous tensor A, is instead symmetric under spin in- 
terchange; it is known to originate from second order 
spin-orbit effects^, but is usually believed to be weak or 
hard to probe. Nevertheless, its electric field dependence 
has not been studied. 

Most interesting effects take place when one of the cou- 
pling coefficients depends linearly on electric field. For 
instance, simple models based on electronic^ or lattice 
mediated polarizatiorP predict that the DM vector is 
electric-field dependent according to D nm oc E x (R„ — 
R m ); this gives rise to the phenomena of magnetically 
induced fcrroclcctricity observed in a large class of ma- 
terials, the cycloidal multiferroics^^. 

In addition, the linear magnetoelectric effect makes 
magnetic excitations electrically dipole active. This gives 





FIG. 1: (Color online) Schematics of low temperature phases 
of 7?Mn03: (a) Collinear sinusoidal phase and (b) Non- 
collinear cycloidal phase. 



rise to the electromagnon, the quasiparticle of the mul- 
tiferroic state^E^. The observation of electromagnons in 
optical experiments provide invaluable clues on the sym- 
metry and magnitude of the magnetoelectric coupling 
present in Eq. ([I]). Moreover, the ability to launch, de- 
tect, and control magnons electrically also hold s prom ise 
for novel applications in information processin g* 19 ! 20 !. 

The observation of magnetically-induced 
ferroelectricity^ and electromagnons^ in the class 
of perovskite manganites i?MnO*3 has made this mate- 
rial the prototype for studies of strong magnetoelectric 
effects. Here R is a rare-earth ion such as Dy, Tb, 
Gd or their mixture, e.g., Gd^Tbi-j;. In the i?Mn03 
family, spins are typically o rder ed with a period incom- 
mensurate with the lattice^Sl. Below the first Neel 
temperature (T = 39 K in DyMn0 3 ), the ground state 
of the Mn spins forms the collinear sinusoidal density 
wave depicted in Fig. [lja). At even lower temperature 
(19 K in DyMnOs), another phase transition takes place 
where the Mn spins order non-collincarly in the cycloid 
ground state shown in Fig.Qb). 

The detection of electromagnons in the cycloidal phase 
of i?MnO*3 led to a surprising observation. Optical ex- 
periments showed that two quite strong electromagnons 
are observed in the cycloidal phase, provided the elec- 
tric field of light was directed along the crystallographic 
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direction This remained true even when the 

cycloid plane was flipped, leading to the conclusion that 
the DM interaction D could not explain the origin of the 
observed strong electromagnon resonances (but in recent 
experiments a weak electromagnon resonance consistent 
with the DM model was observed^) . This is a surprising 
conclusion in view of the fact that the DM interaction is 
known to be responsible for ferroelectricity in these ma- 
terials. 

Optical experiments have also produced a puz- 
zling observation: The lower energy electromagnon, 
unlike the higher energy one, survives also in the 
collinear sinusoidal phase. This is observed, e.g., in 
DyMnOcpS, Gd . 7 Tb .3MnO 3 23 , and Eu^Y.MnOP, 
but in TbMnC>3 no electromagnons are discernible in the 
sinusoidal phased 

Currently, there exists a consensus that the high energy 
electromagnon originat es from magnetostriction, the first 
term in Eq. |l P 3 ! 16 ! 1 !!. However, no consensus exists on 
the origin of the low energy electromagnon. Two quite 
different models were proposed for its explanation: In 
|16j . we showed that magnetostriction plus spin-orbit 
coupling is able to explain the origin of both electro- 
magnons even when the cycloid ground state is purely 
harmonic. In |17) . Mochizuki et al. showed that pure 
magnetostriction plus cycloid anharmonicity (without a 
tensor A) is able to explain the two electromagnons of 
the cycloid phase, suggesting that anharmonicity plays 
a vital role (similar results for BiFe03 were proposed in 
|25|). But neither of the two above-mentioned models is 
able to explain the optical activity of the low energy elec- 
tromagnon in the sinusoidal phase. 



II. MODEL FOR i?Mn0 3 

Here we present a model of electromagnon excita- 
tions that can explain the optical experiments in both 
the sinusoidal and in the cycloidal phases. Our model 
Hamiltonian consists of spin and phonon couplings, H = 



+ H 



P h 



zj(i) , H m 

11 mc i iJ rac ■ 



Here Hs describes exchange 



interactions and single-ion anisotropies, 

H s = 4,mS„ • S m + D a ^(S„ • a) 2 - D b ^(S n • b) 2 . 

n,m n n 

(2) 

We assume D a > and Dt > 0, favoring alignment along 
the b direction. The spins are coupled by exchange inter- 
actions J n ,mi with nearest-neighbor interactions in the 
ab plane denoted by J , next-nearest- neighbor interac- 
tion along b denoted by J2b, and interaction along the 
c direction denoted by J c . The interaction Jo < is 
ferromagnetic while J2b > and J c > are both antifer- 
romagnetic. 

At sufficiently low temperatures, provided that the 
stability condition J^b > —Jo/2 is satisfied, the com- 
petition between the nearest-neighbor ferromagnetic ex- 
change, and the antiferromagnetic next-nearest-neighbor 



exchange favors incommensurate spin ordering. Between 
the first and second Neel temperatures the spins order in 
a sinusoidal density wave, 

5 (R, T) = ±S(T) cos(Q • R + <p)h, (3) 

with S(T) a temperature dependent amplitude [see Eq. 
(3.14) in Ref. [25] for its dependence on model parame- 
ters]. The magnitude of the sinusoidal wave vector Q 
is given by cos(Q6/2) = —J /(2J 2 b). The upper sign in 
Eq. ([3| corresponds to ab layer spins with the integer c- 
coordinate, while the lower sign applies to spins in the 
neighboring ab layers a distance c/2 above and below 
them. 

Our phonon Hamiltonian is 



H p h = 



UJ n X r 



(4) 



where ujq is the (bare) phonon frequency, m* is the effec- 
tive mass, x n is the relative displacement between anions 
and cations in the nth unit cell, e* is the Born charge, 
and E is the electric field of light. 

We divide the linear magnctoclcctric couplings in our 

model into two separate terms, Hml and Hml- The first 
interaction, 
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2,n+a) 



-(1 -> 3,2 -> 4) 



ab 

°l,n+b 



,+a) 



(5) 



does not give rise to electromagnons in a collinear ground 
state, but is necessary to explain the origin of the low 
frequency electromagnon in the cycloidal phased The 
second interaction is instead given by 



Hi 2 } 



X n{[9bc(S\ n 5 l n _|_ 6 )((§2,„ + 52 n +a ) 



9ab(Si n 
(1 ^ 3,2 



+ 4)}, 



2,n 



5*2, n+a) 



(1 o 2)] 
(6) 



where g a b and gb c are coupling constants that can be 
obtained by microscopic calculation (e.g., using density 
functional theory). Like Eq. (IHj), this spin-symmetric in- 
teraction is also invariant under the Pbnm space-group 
operations of i?Mn03, and is therefore consistent with 
lattice symmetry. Both interactions represent anoma- 
lous magnetoelectric coupling, with particular forms of 
the anomalous tensor A [Eq. Q]. A generalization of 
Moriya's theory^ to allow for magnetostriction effects 
shows that such interactions can originate from cross- 
coupling between spin-orbit and magnetostriction effects. 
However, a full microscopic theory is still needed to con- 
firm this expectation. 



III. ELECTROMAGNON SPECTRA 

We adopt the molecular field approximation and ex- 
pand the Hamiltonian H by keeping only terms quadratic 
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in the fluctuation operators, e.g., 5S%, SS a 5P a , 8P%, 
etc. We parametrize the spin excitations 5S = S — Sq 
by 5Si. n = sf n & ± n c, and compute the equations 
of motion using the canonical commutation relations, 
[Sjn>Skni\ ~ i$jk5nmS\- In addition, we also adopted 
the random phase approximation (RPA), i.e., we made 
the substitution S h k -> (S b k ) = S {R k ,T) ■ b in the com- 
mutator above. Such an approximation is expected to 
hold when the fluctuation effects are not too large (i.e., 
we are sufficiently far from the Neel temperature). 

After some manipulation the coupled equations of mo- 
tion for spins and polarization is given by, 



(w 2 - n 2 c J (a% + 




*% + 


S iq) 


= O c , q r 9 , (7a) 


(^-^, q+ko )K- 


4 q + 


b 3q 


Siq) 


= ^Cq+koTg+fco, 










(7b) 


(c 2 - nl c>q ) {s a lq + 


S 2q 


b 3q 


S 4q) 


= 0, (7c) 


- ^EC,q+k ) ( S lq ~ 


S 2q 


4q + 


S 4q) 


= 0, (7d) 



where a = a, a Here sf q is the momentum representation 
of the spin fluctuation sf n . Equations (7a) and (7b) are 
related by a shift in momentum space, q •<-» q + ko where 
fco = 27r/& is the Brillouin zone-edge for magnons. Such 
a relationship corresponds to the fact that "anti-phase" 
fluctuations of neighboring spins with wave vector q are 
equivalent to "in-phase" fluctuations at q + 2ir /b. They 
describe a mode here referred to as a cyclon, with dis- 
persion Qc,q- Similarly, Eqs. (7c) and (7d) share the 
same momentum shift relationship, but describe a differ- 
ent mode referred to as an extra-cyclon. The cyclon and 
extra-cyclon dispersions are shown in Fig. [2] We note 
that the cyclon has a gap proportional to D b , while the 
extra-cyclon has a gap proportional to (2J C + Df,). 

I I I I (2) 

Equations (7a) and (7b) show that H me couples only 
a single electromagnon, the cyclon at q = fco — 2Q, to 
the polar phonon. This takes place through dynamic 
magnetoelectric coupling fJcql^ with 



r„ = 



q 



pafc _ p&c^ pbc _ Y q lb (g ab -> g bc ) y 

g ab S(T) 2 v sin )<5P Q 



(8) 



?q-k a +2Q 



5 9 -fc -2Q + e 2lcl> 5 q+ko 



+2Q 



'2 id 



?q+k Q -2Q\ 




2%/b 



FIG. 2: (Color online) Typical dispersion curves for magnon 
wavevector q along b in the sinusoidal phase: cyclon (black 
solid) and extra-cyclon (red dashed). In the sinusoidal state, 
only the low-energy electromagnon (filled circle) is activated, 
through HSl [Eq. ^]. In the cycloidal state, H$hl [Eq. S] 
activates also the zone-edge electromagnon (hollow square)! 



cyclon at q = fco — 2Q is actually an electromagnon, with 
spectral weight given by 



47TXqCJ 2 A 2 

, ,2 _ q2 

^0 "C,fc -2Q 



(11) 



Here xo — e* 2 /(m*«oWo) is the zero-frequency suscepti- 
bility, with w the unit cell volume. The frequency shift 
A is calculated to be 



A 



S(T)V 2 ( 5bc - 5afc ) 2 tan 2 (f) fi* )feo 



-2Q 



2m* (w§ 



^C,k -2Q> 



) h 4 (?) 



COS H 



J2b 

(12) 

apart from smaller terms of order (g ab — gb c )9ab- Since 
the magnitude of the frequency shift A is the same for 
the magnon and the phonon, we confirm the oscillator 
strength sum rule S em + S p h — 4^X0^0- 



(9) 



Optical experiments such as transmissivity or reflec- 
tivity probe the frequency dependence of the dielectric 
function e(w). After a linear response calculation we ob- 
tain 



IV. ADDITIONAL CONSEQUENCE OF 
ANOMALOUS MAGNETOELECTRIC 
INTERACTION: INCOMMENSURATE 
OSCILLATORY POLARIZATION 



"C,fc - 



2Q 



- A 2 - i 



A 2 -i 



(10) 



Hence e(w) can be written as two Lorentzians, with poles 
at downshifted magnon and upshifted phonon frequen- 
cies. The pole at the magnon frequency shows that the 



In addition to the sinusoidal electromagnon, the cou- 
plings described by Eqs. ([5| and Q have an important 
observational consequence: They lead to an incommensu- 
rate oscillatory polarization (IOP) with wavevector 2Cp^i. 
Minimizing Eqs. |4|-([6]) with respect to the polar phonon 
displacement x„ and plugging in the cycloidal spin order 
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Electromagnons (Far-IR) 


Atomic disp. (X-ray) 


Cycloidal 


(gb + gc) 
(gb - g c ) 


(gb - gc) 
gbc 


Sinusoidal 


(gbc — g a b)g a b 


gb 



TABLE I: This table relates our predictions for optical and 
X-ray experiments to the anomalous magnetoelectric coupling 
parameters g a and g a p introduced in this paper (a, f3 = a,b,c 
are crystallographic directions). The first column refers to 
measurements of electromagnon spectral weight using far-IR 
optical experiments, and the second column refers to the mea- 
surement of magnetically induced lattice distortions using X- 
ray spectroscopy. Each experiment has different magneto- 
electric signatures, depending on whether the ground state 
is cycloidal or sinusoidal. All parameters can be measured 
individually, and in addition, parameters gb and g c can be 
cross-checked. 



we get 

^ = 4 Xo S 2 sm( ( f){(g b ~g c )sm[Qb(2n + l)} 

- gbc cos [Qb(2n + 1)] + g bc } k. (13) 

Note how gbc generates a combination of static and os- 
cillatory polarization along the a direction. 

When the system goes into the sinusoidal phase, this 
polarization changes discontinously to 

^ = 4 X0 S(T) 2 sin g b sin [(2n + l)Qb]k. (14) 

Such an oscillatory polarization can be detected by X-ray 
scattering. Indeed, Kimura et at detected an oxygen 
oscillation with wavevector 2Q in both the cycloidal and 
sinusoidal phases [see blue dots in Fig. 1(c) of [5j]- Just 
like our prediction, the X-ray intensity in [S] showed an 
apparent discontinuity in oxygen displacements in the 
transition from cycloidal to sinusoidal phase. 

Table I shows how a combination of optical and 
X-ray scattering experiments are capable of measur- 
ing the magnetoelectric coupling constants individually, 
and even cross-check some of them. For instance, for 
DyMnC>3, we obtain from the measured electromagnon 
spectral weights^ and X-ray diffraction intensities^ 
in the sinusoidal and cycloidal phases the values of 
g b ~ 170 erg/ (cm esu), g c ~ -40 erg/ (cm esu), g ab ~ 



100 erg/(cm esu) 3> gbc- In TbMnC>3, the sinusoidal elec- 
tromagnon could not be observed experimentally^. This 

(2) 

indicates that g a bi9bc <S .96, g c , i-e., i?mc is much weaker 
than H^l in TbMn0 3 . 



V. DISCUSSION AND CONCLUSION 

We now consider the justification of our model and 
other possibilities for the activation of the electromagnon 
in the collinear sinusoidal phase. First, we note that the 
presence of DM interaction in principle also predi cts an 
electromagnon in the collinear sinusoidal stat cj 28 * 2 ^ . How- 
ever, this scenario is ruled out by the experiments where 
the electromagnon is activated with E along a only. In 
addition to the DM interaction, there are no other lin- 
ear magnetoelectric couplings anti-symmetric under the 
exchange of spins that would be allowed by lattice sym- 
metry 

The anharmonic cycloid model of Mochizuki et a/P^ 
would give rise to no electromagnon activity in the sinu- 
soidal phase. One possibility for the activation of sinu- 
soidal electromagnons in this scenario would be to include 
additional single-ion anisotropy so that the spins in the 
ground state are tilted off the b axis®]. However, in this 
case both high and low energy electromagnons get acti- 
vated. We found no scenario where deformation of the 
sinusoidal ground state activates the low-energy electro- 
magnon without activating the high-energy one. 

Concerning other possible symmetry-allowed magne- 
toelectric interactions, we note that the other couplings 
quadratic in spin do not couple electric field linearly to 
magnons. More specifically, terms of the form x^SfSj, 
x^SfSj, and x^SfSj lead to contributions that are third 
order in fluctuation operators, and the term x^S^S*- 
does not couple polarization to magnons. Hence in the 
collinear sinusoidal state only the couplings linear in S b 
considered in the present work [Eq. ([6])] can be responsi- 
ble for the electromagnons in collinear sinusoidal state. 

In conclusion, we showed that anomalous magnetoelec- 
tric coupling gives a natural explanation for the origin 
of electromagnons in both the cycloidal and sinusoidal 
phases of i?Mn03. It remains an open question to study, 
e.g., through ab-initio methods, the microscopic mecha- 
nism of anomalous magnetoelectric coupling. 
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